Effective variant filtering and expected candidate variant yield in studies of rare human disease.

In studies of families with rare disease, it is common to screen for de novo mutations, as well as recessive or dominant variants that explain the phenotype. However, the filtering strategies and software used to prioritize high-confidence variants vary from study to study. In an effort to establish recommendations for rare disease research, we explore effective guidelines for variant (SNP and INDEL) filtering and report the expected number of candidates for de novo dominant, recessive, and autosomal dominant modes of inheritance. We derived these guidelines using two large family-based cohorts that underwent whole-genome sequencing, as well as two family cohorts with whole-exome sequencing. The filters are applied to common attributes, including genotype-quality, sequencing depth, allele balance, and population allele frequency. The resulting guidelines yield ~10 candidate SNP and INDEL variants per exome, and 18 per genome for recessive and de novo dominant modes of inheritance, with substantially more candidates for autosomal dominant inheritance. For family-based, whole-genome sequencing studies, this number includes an average of three de novo, ten compound heterozygous, one autosomal recessive, four X-linked variants, and roughly 100 candidate variants following autosomal dominant inheritance. The slivar software we developed to establish and rapidly apply these filters to VCF files is available at https://github.com/brentp/slivar under an MIT license, and includes documentation and recommendations for best practices for rare disease analysis.

CaBagE: A Cas9-based Background Elimination strategy for targeted, long-read DNA sequencing.

A substantial fraction of the human genome is difficult to interrogate with short-read DNA sequencing technologies due to paralogy, complex haplotype structures, or tandem repeats. Long-read sequencing technologies, such as Oxford Nanopore's MinION, enable direct measurement of complex loci without introducing many of the biases inherent to short-read methods, though they suffer from relatively lower throughput. This limitation has motivated recent efforts to develop amplification-free strategies to target and enrich loci of interest for subsequent sequencing with long reads. Here, we present CaBagE, a method for target enrichment that is efficient and useful for sequencing large, structurally complex targets. The CaBagE method leverages the stable binding of Cas9 to its DNA target to protect desired fragments from digestion with exonuclease. Enriched DNA fragments are then sequenced with Oxford Nanopore's MinION long-read sequencing technology. Enrichment with CaBagE resulted in a median of 116X coverage (range 39-416) of target loci when tested on five genomic targets ranging from 4-20kb in length using healthy donor DNA. Four cancer gene targets were enriched in a single reaction and multiplexed on a single MinION flow cell. We further demonstrate the utility of CaBagE in two ALS patients with C9orf72 short tandem repeat expansions to produce genotype estimates commensurate with genotypes derived from repeat-primed PCR for each individual. With CaBagE there is a physical enrichment of on-target DNA in a given sample prior to sequencing. This feature allows adaptability across sequencing platforms and potential use as an enrichment strategy for applications beyond sequencing. CaBagE is a rapid enrichment method that can illuminate regions of the 'hidden genome' underlying human disease.

To ERV Is Human: A Phenotype-Wide Scan Linking Polymorphic Human Endogenous Retrovirus-K Insertions to Complex Phenotypes.

Approximately 8% of the human genome is comprised of endogenous retroviral insertions (ERVs) originating from historic retroviral integration into germ cells. The function of ERVs as regulators of gene expression is well established. Less well studied are insertional polymorphisms of ERVs and their contribution to the heritability of complex phenotypes. The most recent integration of ERV, HERV-K, is expressed in a range of complex human conditions from cancer to neurologic diseases. Using an in-house computational pipeline and whole-genome sequencing data from the diverse 1,000 Genomes Phase 3 population (n = 2,504), we identified 46 polymorphic HERV-K insertions that are tagged by adjacent single nucleotide polymorphisms (SNPs). To test the potential role of polymorphic HERV-K in the heritability of complex diseases, existing databases were queried for enrichment of established relationships between the HERV-K insertion-associated SNPs (hiSNPs), and tissue specific gene expression and disease phenotypes. Overall, hiSNPs for the 46 polymorphic HERV-K sites were statistically enriched (p < 1.0E-16) for eQTLs across 44 human tissues. Fifteen of the 46 HERV-K insertions had hiSNPs annotated in the EMBL-EBI GWAS Catalog and cumulatively associated with >100 phenotypes. Experimental factor ontology enrichment analysis suggests that polymorphic HERV-K specifically contribute to neurologic and immunologic disease phenotypes, including traits related to intra cranial volume (FDR 2.00E-09), Parkinson's disease (FDR 1.80E-09), and autoimmune diseases (FDR 1.80E-09). These results provide strong candidates for context-specific study of polymorphic HERV-K insertions in disease-related traits, serving as a roadmap for future studies of the heritability of complex disease.

Allergies and Childhood Acute Lymphoblastic Leukemia: A Case-Control Study and Meta-analysis.

Background: Allergic disease is suspected to play a role in the development of childhood acute lymphoblastic leukemia (ALL). Studies conducted over the last several decades have yielded mixed results.Methods: We examined the association between allergy, a common immune-mediated disorder, and ALL in the California Childhood Leukemia Study (CCLS), a case-control study of 977 children diagnosed with ALL and 1,037 matched controls (1995-2015). History of allergies in the first year of life was obtained from interviews, mainly reported by mothers. Logistic regression analyses were conducted to estimate ORs and 95% confidence intervals (CIs), controlling for birth order, daycare attendance, and mode of delivery. In addition, we conducted meta-analyses with data from the CCLS and 12 published studies and employed a new method to estimate between-study heterogeneity (R_b).Results: Overall, no associations were observed between childhood ALL risk and specific allergy phenotypes or any allergy, as a group. However, having any allergy was associated with an increased risk of ALL among the youngest study participants. In the meta-analysis random-effects models, reduced odds of ALL were associated with hay fever (metaOR = 0.65; 95% CI, 0.47-0.90); however, restricting the analysis to studies that used medical records for assessment of allergy or recently published studies led to null or attenuated results.Conclusions: Overall, our findings do not support a clear association between allergy and childhood ALL.Impact: The degree to which epidemiologic studies can inform the relationship between allergies and risk of childhood ALL is limited by R_b. Cancer Epidemiol Biomarkers Prev; 27(10); 1142-50. ©2018 AACR.

In utero cytomegalovirus infection and development of childhood acute lymphoblastic leukemia.

It is widely suspected, yet controversial, that infection plays an etiologic role in the development of acute lymphoblastic leukemia (ALL), the most common childhood cancer and a disease with a confirmed prenatal origin in most cases. We investigated infections at diagnosis and then assessed the timing of infection at birth in children with ALL and age, gender, and ethnicity matched controls to identify potential causal initiating infections. Comprehensive untargeted virome and bacterial analyses of pretreatment bone marrow specimens (n = 127 ALL in comparison with 38 acute myeloid leukemia cases in a comparison group) revealed prevalent cytomegalovirus (CMV) infection at diagnosis in childhood ALL, demonstrating active viral transcription in leukemia blasts as well as intact virions in serum. Screening of newborn blood samples revealed a significantly higher prevalence of in utero CMV infection in ALL cases (n = 268) than healthy controls (n = 270) (odds ratio [OR], 3.71, confidence interval [CI], 1.56-7.92, P = .0016). Risk was more pronounced in Hispanics (OR=5.90, CI=1.89-25.96) than in non-Hispanic whites (OR=2.10 CI= 0.69-7.13). This is the first study to suggest that congenital CMV infection is a risk factor for childhood ALL and is more prominent in Hispanic children. Further investigation of CMV as an etiologic agent for ALL is warranted.

Genotyping Oral Commensal Bacteria to Predict Social Contact and Structure.

Social network structure is a fundamental determinant of human health, from infectious to chronic diseases. However, quantitative and unbiased approaches to measuring social network structure are lacking. We hypothesized that genetic relatedness of oral commensal bacteria could be used to infer social contact between humans, just as genetic relatedness of pathogens can be used to determine transmission chains of pathogens. We used a traditional, questionnaire survey-based method to characterize the contact network of the School of Public Health at a large research university. We then collected saliva from a subset of individuals to analyze their oral microflora using a modified deep sequencing multilocus sequence typing (MLST) procedure. We examined micro-evolutionary changes in the S. viridans group to uncover transmission patterns reflecting social network structure. We amplified seven housekeeping gene loci from the Streptococcus viridans group, a group of ubiquitous commensal bacteria, and sequenced the PCR products using next-generation sequencing. By comparing the generated S. viridans reads between pairs of individuals, we reconstructed the social network of the sampled individuals and compared it to the network derived from the questionnaire survey-based method. The genetic relatedness significantly (p-value < 0.001) correlated with social distance in the questionnaire-based network, and the reconstructed network closely matched the network derived from the questionnaire survey-based method. Oral commensal bacterial are thus likely transmitted through routine physical contact or shared environment. Their genetic relatedness can be used to represent a combination of social contact and shared physical space, therefore reconstructing networks of contact. This study provides the first step in developing a method to measure direct social contact based on commensal organism genotyping, potentially capable of unmasking hidden social networks that contribute to pathogen transmission.

Mode of delivery and risk of childhood leukemia.

BACKGROUND: Childhood infection and immune response have long been suspected in the etiology of childhood leukemia, specifically acute lymphoblastic leukemia (ALL). Normal primary inoculation of the core human microbiome is circumvented by cesarean section (CS) delivery, which is a proposed modulator of both immune response and early-life infection. METHODS: In this study, we examined CS delivery and the risk of childhood leukemia using data from the California Childhood Leukemia Study (CCLS) case-control study and additive logistic regression models. RESULTS: We observed no association between CS and acute myelogenous leukemia [OR, 0.96; 95% confidence interval (CI), 0.52-1.55]. We observed a suggestive association for ALL and CS (OR, 1.22; 95% CI, 0.97-1.54). When examining common ALL (cALL), defined as ALL with expression of CD10 and CD19 surface antigens and diagnosis occurring between 2 and 5.9 years of age, we found a significant association with CS (OR, 1.44; 95% CI, 1.0-2.06). ALL subjects that are not cALL showed a similar risk as ALL overall (OR, 1.15; 95% CI, 0.91-1.44). Because of previous findings suggesting effect modification, we stratified cALL subjects by Hispanic status. Although we observed no relationship for CS in non-Hispanics (OR, 1.14; 95% CI, 0.72-1.79), we did observe a strong association between cALL and CS in Hispanics (OR, 2.34; 95% CI, 1.23-4.46). CONCLUSION: Within the CCLS, CS delivery seems to be associated with cALL and Hispanic subjects may be driving the association. IMPACT: Further research combined with investigations into response to early infection and the microbiome is warranted.